function ret=my_sqrt(A)
 vect=getMainVector(A);
 vect=vect./2;
 mat=u_rot(vect);
 s=-mat(1,2)*mat(2,1);
 theta=asin(sqrt(s));
 theta=theta*2;
 if(s==1)
  c=(mat(1,2)-mat(2,1))/(2*sin(theta/2));
  sub=acos(c);
  ret=[sub*2 theta 0];
  return;
 endif
 if(theta==0)
  error("Zero divide exception. Other approach will solve");
 endif
 c=(mat(1,2)-mat(2,1))/(2*sin(theta/2));
 sub=acos(c);
 c=(mat(1,1)+mat(2,2))/(2*cos(theta/2));
 sum=acos(c);
 ret=[sum+sub theta sum-sub];
endfunction


